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Abstract. 

We present a general model of disorder in Kondo alloys that, under certain 
conditions, leads to non-Fermi liquid behavior. The central underlying idea is the 
presence of a distribution of local Kondo temperature scales. If this distribution is 
broad enough, such that there are sites with arbitrarily low Kondo temperatures, a 
non-Fermi liquid phase is formed. We analyze thermodynamics and transport in this 
approach and show it is consistent with a number of Kondo alloys. We also compare 
the predictions of this model with the measured dynamical magnetic response of 
these systems. 



1. Introduction 

Since its introduction by Landau in 1956, the Fermi liquid paradigm has been the 
most robust stepping-stone of the theory of metals [Q. Its central assumption is the 
existence of a one to one correspondence between the excitations of a free Fermi gas 
and an interacting fermionic system. Even though originally formulated to explain 
the behavior of a neutral Fermi system ('^He), it has formed the basis of much of 
our understanding of what goes on in metallic systems. It lies behind such successful 
theories as Migdal's electron-phonon theory and provides the background upon which 
instabilities such as superconductivity can be understood. It has proved to be a 
valuable guide in the analysis of even very strongly correlated systems, such as some 
of the heavy fermion metals, where the theory seems to survive the most extreme 
circumstances, with enormous renormalizations in the quasiparticle effective mass. 

The very idea of a correspondence between the low-lying excitations of an 
interacting system and those of a reference one has proved more general than the 
original strict application envisioned by Landau. Indeed, this principle of "adiabatic 
continuity" [|j can be found in very diverse situations with very different and sometimes 
non-trivial reference systems. For example, in the so-called Fermi liquid theory of 
the Kondo or Anderson impurity problem, the reference system is a non-interacting 
scattering center embedded in a free Fermi sea[^ Q . Other examples include superfluid 
■^Hel^, atomic nuclei||^ and interacting disordered metalsQ]. The corresponding 
reference systems, in these latter cases, are the BCS model of pairing, a shell model 
of the nucleus and a model of diffusive electrons in a random potential, respectively. 
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Given this operational definition of a Fermi liquid, the question of whether a 
particular system behaves as a Fermi liquid or not can be a rather delicate one. It 
is necessary first to determine what is a natural reference system for a particular 
compound and then find out whether the compound can fit a description which is 
adiabatically connected to the reference one. We cannot overemphasize the fact that 
this is not always a straightforward task. Daunting as the task may be, researchers 
have nonetheless developed rules of thumb to classify various systems as non-Fermi 
liquids. One such commonly used criterion is a resistivity which depends linearly on 
the temperature at low temperatures, a result certainly inconsistent with a clean Fermi 
liquid. However, as helpful as these empirical tests are, as a matter of scientific rigor, 
we must stress that the classification of a certain system as a non-Fermi liquid often 
involves careful cross-checking between experiment and theory. This is particularly 
critical when the system is disordered. With these caveats in mind, we can proceed 
to consider the question of non-Fermi liquid behavior. 

In recent years, a growing number of metallic compounds have come to be known as 
counter-examples to the old Fermi liquid paradigm. This is usually taken to mean that 
the adiabatic continuity hypothesis appears to break down, if one assumes the most 
natural reference system for the system under study. This has been perhaps most 
strongly emphasized in the case of the high Tc cuprate superconductors!^, where 
Landau's original Fermi liquid theory is unexpectedly violated in the normal state. In 
addition, some heavy fermion systems have also been discovered which do not seem 
to fit the general picture of a Fermi liquid |lOj. The bulk of this evidence has in 
effect turned the study of non-Fermi liquid behavior into an independent frontier area 
of research whose ramifications can only be dimly glimpsed. The main goal is to 
determine the possible routes towards the breakdown of Fermi liquid theory and to 
classify what new low-lying excitations are present in these systems. 

Among the non-Fermi liquid candidates, some should be classified as belonging to 
the heavy fermion family. This is taken to mean, in this context, that for some range of 
composition or for certain values of external parameters - such as pressure or magnetic 
field - these compounds show typical heavy fermion physics with high temperature 
incoherent Kondo behavior accompanied by the formation of a low temperature heavy 
Fermi liquid, with or without magnetic or superconducting order[^. However, for 
the more interesting values of these parameters, their behavior is not that of a Fermi 
liquid. 

Some of the heavy fermion non-Fermi liquid metals have been convincingly 
associated with the proximity to a quantum critical point. These are the 
alloy CeCu6-a;Au2:||l^ and the compound CePd2Si2||l^. Both systems show 
antiferromagnetism: CeCug-xAux for a; > 0.1 and CePd2Si2 for pressures P < 26 
kbar. At the critical values of these parameters, the Neel temperature T/v vanishes 
and the behavior of the system appears to be governed by the zero temperature critical 
point. The anomalous non- Fermi liquid behavior of these compounds is exemplified by 
the resistivity which is given by p « po + AT in CeCus.gAuo.i and as p « po + BT^'^ 
in CePd2Si2 at P = 26 kbar. In the latter case, the non-Fermi liquid behavior is 
interrupted by what appears to be a phase transition into a superconducting state at 
Tc « QAK . The thermodynamic response of CeCus.gAuo.i is also anomalous with the 
specific heat C/T w aln(r/ro) and x ~ Xo(l — aVT). In addition, external pressure 
can suppress T/v to zero in CeCus.rAuo.s and the anomalous behavior associated with 
the quantum critical point can be recovered. Though the non-Fcrmi liquid behavior 
can be reasonably ascribed to the proximity of a quantum critical point a complete 
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theory does not yet seem to exist [Q. 

In contrast to the latter compounds, a series of other alloys also seems to 
show characteristic non-Fermi liquid behavior which cannot be clearly associated 
with quantum criticality. A partial list of these is given in Table ^ together with 
their corresponding resistivity, specific heat and magnetic susceptibility (see also 
Refs. AH of them show anomalous thermodynamic and transport properties 

incompatible with a Fermi liquid description. 



Table 1. Heavy fermion alloys which exhibit non- Fermi liquid behavior not obviously 
ascribed to the proximity to a quantum critical point and their properties. Here, p{T) 
is the DC resistivity, C{T) is the specific heat, xC^) the magnetic susceptibility 
and 1/t{u!) is the frequency dependent scattering rate. Below, A> and luq > 0. 



Compounds 



PiT) 



C{T)/T x{T) 



l/r{u)(T = Q) 



UCu5-.Pd,ll|,|l§ 
Mi_,U.Pd3 (M=Sc,Y)0 @ 
Lai_3:Ce3;Cu2.2Si2 ||li| 
Ui-^Th.PdaAlajTQ 
Cei-^j^Th^RhSbpl 



Po-AT 
po - AT 
Po - AT 
Po-AT 



aln{To/T) a\n(To/T) 

a\n{To/T) AT-°''^ 

a\n{To/T) aln{To/T) 

aln(To/r) xo - AVT 

a\n{To/T) 



l/ro(l ■ 
l/^o(l 



- (jj/uJq) 



U.Thi-.MaSia (M=Ru,Pt,Pd)[M po + ^InT aln(ro/r) aln(ro/r) 



There are a few proposed scenarios to try to explain this anomalous behavior. 
Some rely on the existence of a critical point at T = OQ. Though these theories 
should be relevant to the cases of CeCue-icAua: and CePd2Si2 alluded to above, the 
alloys in Table |l] are not obviously close to a phase boundary. 

On the other hand, other proposals have focused on a local approach to the non- 
Fermi liquid physics. These include exotic impurity models governed by a non- Fermi 
liquid fixed point such as the quadrupolar Kondo model and various multi-channel 
Kondo models!^, These are largely based on the anomalous behavior of a dilute 
system of such impurities. The inclusion of lattice effects presents an additional 
challenge that has only recently started to be addressed Another proposed route 
to non- Fermi liquid behavior has been to consider the competition between local charge 
and spin fiuctuations[^. It is quite possible, perhaps even likely, that the origin of the 
anomalous behavior is different for different systems. We mention, in particular, the 
case of the last entry of Table |], U2;Thi_i,M2Si2 (M=Ru,Pt,Pd), whose description 
as a dilute system of magnetic two-channel impurities |^ in Ref. appears to be 
fairly good. 

The fact that the systems in Table ^ are all disordered alloys immediately poses 
the question of the role of disorder in the formation of the non-Fermi liquid state. 
Furthermore, recalling our operational definition of a Fermi liquid, one must not 
neglect the fact that the reference system in this case is almost certainly a disordered 
one. The behavior of local moments in disordered systems has been considered in past 
studies H, ||. 

In this paper, we will focus on the first alloy of the table, UCu5_a;Pd2;. The 
accumulated experimental data on this system has suggested to us that a model 
of disorder in f-electron systems is able to account for its anomalous behavior. 
Paramount to this conclusion was the Cu NMR study of Ref. [Q, which detected 
a large inhomogeneous broadening of the NMR line, attributable to microscopic 
disorder. It is the goal of this paper to present a scenario in which a non-Fermi 
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liquid state is generated as a consequence of the interplay of disorder and strong 
correlations . This theory provides a consistent way of understanding the non- Fermi 
liquid anomalies in both transport and thermodynamic properties based on a single 
underlying mechanism. The central idea of this theory is that moderate bare disorder 
in a lattice model of localized moments is magnified due to the strong local correlations 
between the moments and the conduction electrons. In particular, a broad distribution 



of local energy scales (Kondo temperatures) is generated! 29 . A few local sites with 



very small Kondo temperatures are unquenched at low temperatures and dominate the 
thermodynamics and transport, forming a dilute gas of low- lying excitations above the 
disordered metallic ground state. The presence of these unquenched moments leads to 
the formation of a non-Fermi liquid phase. While we think that VCu^-xP^ix (a; = 1 
and X = 1.5) provides the best candidate system to exhibit such a phase to this date, 
we expect this type of behavior to be seen in other Kondo alloys. Whether the other 
alloys in Table |l] can be understood in this framework is not clear at this moment. 

The outline of this paper is as follows. In Section H, we discuss the experimental 
work on the VGn^-xPdx alloys which served as the motivation for this study, 
emphasizing the most important underlying physical ideas. In Section H, we present 
the model and our dynamical mean field theory approach to its solution. In Section ^, 
we show the nature of the ground state of the model in the clean as well as in the 
disordered limits. In Section we explain the origin of the linear dependence on 
the temperature of the resistivity. In Section ^ we apply this disorder model to the 
description of the dynamical susceptibility and compare the results with experiments 
on UCu5_a;Pda;. Finally, we conclude with a discussion of the limitations of our 
approach in Section ^ In order to lighten the line of argument, a few derivations are 
left to the appendices. 



2. Experimental basis for the disorder model 

The major reason for considering disorder as the possible origin of the anomalous 
behavior of the alloys of Table |l| was the Cu NMR study reported in Ref. ||30|| . 
In this paper, the Cu NMR field-swept powder pattern spectra of UCu5_a;Pda: for 
X = I and X = 1.5 were shown to exhibit strong inhomogeneous broadening that 
could only be explained by invoking the presence of short-range disorder. A simple 
disorder model was then used to describe both the broadening of the NMR line and the 
spin susceptibility and specific heat measured in these samples. The model consisted 
of a collection of independent spins, mimicking the Uranium ions, each coupled to 
the conduction electron bath by a dimensionless Kondo coupling constant A = poJ, 
which was allowed to be randomly distributed in the samples. These were supposed 
to mimic the local disorder induced by the Pd substitution in the Cu ligand sites 
of the parent UCu5 compound and, for simplicity, were assumed to be distributed 
according to a Gaussian. Assuming that the spins were completely uncorrelated, the 
thermodynamic response was then calculated by taking an average over the response of 
a single Kondo spin with the distribution of coupling constants. Since the physics of a 
single Kondo spin is characterized by a single energy scale Q - the Kondo temperature 
- the important quantity in the averaging procedure is the distribution of Kondo 
temperatures. 

The Kondo temperature is defined as 



Tk = De-^/^. 



(2.1) 
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Because of the exponential dependence on the couphng constant A, the corresponding 
distribution of Kondo temperatures is skewed and considerably broadened. Figure ^ 
shows the experimentally determined distributions of Kondo temperatures for both 
alloys (a; = 1 and x ~ 1.5). We first note that, due to the Jacobian in the definition 
of P{Tk), it actually diverges weakly as Tk 

P(T.) . P[A(T.)] (^) = ^ r^l^'i... . (2.2) 



Tk (HD/Tk)) 




I 100 200 300 400 
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Figure 1. Experimentally determined distributions of Kondo temperat ure s of the 

alloys UCu5_:iPda: with x = 1( ) and 1.5( ) (from Ref. [Q). The 

shaded area below T represents the low-T^ spins which remain unquenched at that 
temperature. The upturn at very low Tx's is not shown (see text for discussion). 

However, for the distributions determined experimentally for UCus-^^Pda; (a; = 1 
and X = 1.5) and shown in Figure ^ the point at which P{Tk) shows an upturn 
and starts to diverge occurs at very low temperatures: Tk ~ 2AK {x — 1) and 
Tk ~ 0.8A' {x — 1.5), which can hardly be distinguished on the scale of Figure |l|. 
None of the conclusions derived from the use of the full P{Tk) below depends on this 
small diverging tail, as long as one does not probe into the low temperature region 
below the upturn point scale. Instead, as will be shown below, the response of the 
system is dominated by the part of P{Tk) above the upturn, which appears to be 
tending to a constant value as Tk — > 0. This means that there is a finite number of 
spins with arbitrarily low Tr-'s in the sample, a feature that will be at the root of 
the non-Fermi liquid features. We have, therefore, chosen to depict P{Tk) only above 
the upturn point and will confine the analysis below to a distribution which tends 
to a constant as Tk — > 0. It remains an interesting possibility whether the effects of 
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the low Tk divergence of P{Tx) can be actually observed. On the other hand, it is 
also conceivable that there is a physical infrared cutoff Ac to P(A) (Ac ~ 0.11 for both 
distributions) . In the absence of any experimental evidence for either possibility, we 
will focus on the more relevant effect of the distribution of Kondo temperatures above 
the upturn. 

From the solution of the single impurity Kondo problem it is known that 
thermodynamic quantities like the impurity spin susceptibility x{T) increase with 
decreasing temperature in a Curie fashion, with logarithmic corrections However, 
at very low temperatures, the Curie-like divergence is cut off and the susceptibility 
saturates to a constant value proportional to the inverse Kondo temperature. The 
scale that separates the high temperature from the low temperature region is the 
Kondo temperature. The impurity specific heat divided by the temperature Cy(T)/T 
has a somewhat similar behavior. The saturation at the lowest temperatures is a 
consequence of the "disappearance" of the free spin response, through the formation 
of a singlet ground state with the conduction electron bath, a process generally known 
as "quenching". 

If we are measuring the thermodynamic properties at a certain finite temperature 
T, there will always be Uranium ions with Tk < T, which remain unquenched and 
whose contribution dominates the overall response (the shaded area in Figure |l|) . As 
the temperature is lowered, the number of such spins decreases, as more and more of 
them become quenched. The thermodynamic behavior of the disordered system is, 
therefore, dominated by the tails of the distribution of Kondo temperatures, rather 
than by the average, a situation commonly known as a Griffiths phase . 

Since both x{T) and Cv{T)/T scale as the inverse Kondo temperature at T = 0, 
the fact that P{Tk = 0) 7^ immediately implies that the leading behavior of the 
averaged quantity is a logarithmic divergence. Indeed, consider, for example, the 
susceptibility 

where the asymptotic forms of f{x) are known 

{a - Px^ X < 1; 

X \ mx J 

a, (3 and 7 being universal numbers. To find the leading low temperature behavior it 

is sufficient to use the first term in P{Tk) = Pq + PiTr H • Therefore, if (• • Y'^ 

denotes the average over the distribution of Kondo temperatures 

{yc{T)r ^ ^P{TK)f 

^P,f{l/y), (2.5) 

y 

where we cut off the integral by an arbitrary scale F which sets the region of validity 
of the approximate form for P{Tk)- From ( ^.4|) , it is clear that the lower limit of the 
integral in (2.5) gives a regular contribution whereas the upper limit dominates 

{x{T)r « -aPo ^ «Poln (I) . (2.6) 



y \T 
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A similar analysis holds for Cv{T)/T, with a similar divergence. 

A caveat about the experimental situation is in order here. In the literature, one 
often finds different power laws fitted to thermodynamic and transport properties 
of this and other compounds. In the particular case of UCu5_a;Pd2:, one can find 
X(T) = xoT-'', with 77 = 0.27 ± 0.03|l|, = 0.25[| and 77 = I/3[|§ ||. It 
is not always clear what temperature range was used in the fits. In the case of 



Ref. 1 24 , the t] = 1/3 power law is argued to be valid in an intermediate range of 
temperatures (between 20 and 300K). Besides, there is a small sample dependence to 
these quantities and, in the case of the susceptibility, the magnetic field strength used 
in the measurement can have an important effect We feel that all these aspects 
should be carefully considered when trying to determine the temperature dependence. 

It is important to notice that the distribution of Kondo temperatures reported 
in Ref. which is nearly featureless at low Tk, gives a logarithmic divergence 

in xiT) and Cv(T)/T as the leading low temperature behavior. In the fits shown 
in Ref. ||30| , there are clear deviations from this leading behavior at intermediate and 
higher temperatures. We point out that power law behavior could be inferred from the 
analysis of a narrow window of temperatures. What gives us particular confidence in 
the disorder model is the fact that a full theoretical curve, with the added complication 
of a finite magnetic field, can be well fitted to the experiments. On the other hand, it 
is not clear how changes in the specific form of the distribution function used would 
affect the quality of the fit in the intermediate and high temperature regions. 

Finally, for completeness, we have plotted in Figure ^ the Wilson ratio 
{Rw = Tx(T)/CviT)) prediction in the disorder model as a function of temperature. 
When P(0) = 0, the Wilson ratio tends to its universal value of 2 when T ^ 0, which is 
consistent with the Fermi liquid prediction. When P{0) 7^ 0, the Wilson ratio appears 
to tend to a different value for a wide range of temperatures. This can be understood 
from Equation ( ^.6| ). Because the scale F, which enters this equation, depends on 
the details of the full scaling curve, it is going to be, in general, different for xC^) 
compared to Cv{T)/T. Though the scale is asymptotically irrelevant when T ^ F, 
because of the slow logarithmic dependence, it takes a temperature that is too low in 
order for one to be able to observe the asymptotic behavior. 

We point out that one feature which distinguishes the alloy UCus-a^Pda, from the 
other alloys in Table [| is the nature of the doping. In the former case, the f-sublattice 
remains unchanged and substitutions are introduced in the ligand sites. In the other 
cases, the doping is performed directly in the f-sites. It is apparent that ligand site 
substitution will affect the hybridization matrix element V of an Anderson model 
description and that will likely be the primary effect (see Section ^ for a discussion 
of the disordered Anderson lattice model). It is less certain what kind of effect such 
substitution will have on the f-site energy Ef. Moreover, lattice microstrains should 
also be reflected in the hybridization amplitude V. This is probably what happens in 
UCu4Pd, which could be stoichiometric and indeed appears to be so according the X- 
rays, but whose broad NMR lines give unequivocal evidence of disorder. On the other 
hand, f-site replacements can be mimicked by a different f-site energy Ef {Ef ^ +00, 
when the substitutional ion is a non-magnetic "Kondo hole" ) and, presumably, a 
different hybridization V. In the interest of generality, and in the absence of a better 
understanding of the nature of the disorder, we will consider in this paper both Ef 
and V disorder. What should be emphasized, however, is the fact that, whatever the 
microscopic origin of disorder, its observable consequences will always depend only on 
the distribution of Kondo temperatures P{Tk), for it is these scales which govern the 
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Figure 2. The Wilson ratio as a function of temperature as predicted by the disorder 

model for the cases where P{Tk = 0) ^ ( ) and P{Tk = 0) = { ). Note 

how the latter tends to the universal zero temperature value of 2. 



measured responses. In other words, almost all the information that can be obtained 
experimentally about the disorder will be filtered through P{Tk)- 

Therefore, the main conclusion of Rcf. |30| was that a simple model of Kondo 
disorder is able to account for the NMR line broadening as well as the anomalous 
thermodynamic response of the alloys studied. This simple model relied on the 
assumption that the correlations between the Kondo spins could be neglected, though 
they form a concentrated lattice in these alloys. Past studies of heavy fermion 
systems, both experimental and theoretical! 35 , give us confidence that this may 



not be a bad approximation for thermodynamic quantities (see also the Appendix C). 
However, it is well known that transport is very different in dilute Kondo systems as 
opposed to concentrated Kondo lattices. The so-called onset of coherence in the clean, 
concentrated case is the most dramatic example. It is characterized by rather large 
resistivities at high and intermediate temperatures, consequence of strong incoherent 
Kondo scattering off the localized moments, which is then followed by a precipitous fall 
by some orders of magnitude at the lowest temperatures. At very low temperatures, 
the system can then be characterized as a good metal. This low-temperature coherence 
is a consequence of translational invariance and Bloch's theorem in the ordered lattice 
system. We will, therefore, address the question of whether a single unifying approach 
to the disorder problem, sensitive enough to the formation of a coherent state in 
the clean limit, is able to be formulated. We will answer the latter question in the 
affirmative and will thereby show that the same mechanism that leads to diverging 
thermodynamic properties at low T - the presence of spins with arbitrarily low T^s 
- will also predict a resistivity which is linear in T, as observed in many non-Fermi 
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liquid heavy fermion alloys. 



3. The disordered Anderson lattice model and the dynamical mean field 
theory equations 

In order to consider the interplay of disorder and local moment behavior we will focus 
our attention on a model of disordered Anderson lattices. We will be guided by the 
series of alloys listed in Table |l|. In the spirit of previous works on heavy fermion 
systems we will assume that the f-electron sites can be described by simplified non- 
degenerate Anderson impurities hybridized with a single broad uncorrelated band of 
conduction electrons 

H = J2 <k)cy~^ + Eff]j^^ + Y.^V,c]j^^ + H.c.) + C/^n/.Tn/,^, (3.1) 

k<j 3<^ 3<^ i 

where, c^^ destroys a conduction electron with momentum k and spin a from a band 
with dispersion e(fc) and half bandwidth D and fja destroys an f-electron at site j 
with spin a. Since U is generally the largest energy scale in typical f-shell parameters 
it will be taken to infinity in this paper. 



The important thing to notice in Equation (3.1) is that, unlike in the usual periodic 
Anderson model, the local f-shell parameters Ej and Vj are taken here to be random 
numbers distributed, in general, according to two different distributions Pi{E^) and 
P2{V). As explained before, in the absence of a more detailed understanding of the 
microscopic nature of the disorder in these systems, one should take P\{Ef) and 
P2{y) to be, in principle, of the most general form. Note, however, that we will not 
assume the disorder widths to be too large. For instance, the NMR study of Ref. 
suggests bare distribution widths to be around 20% of their average values. However, 
as we will see, correlation effects will themselves generate large effective disorder. In 
practice, we have considered Gaussian and uniform distributions. Disorder in the 
conduction electron band, though certainly present in these alloys, is not subject to 
these correlation-induced renormalizations and has therefore been neglected in our 



treatment. Despite these simplifications, we believe the model in (3T) captures the 
essential ingredients relevant to the study of disorder in Kondo alloys. 

Some previous studies of disordered Anderson or Kondo lattices have been reported 
in the literature. We mention, in particular, the early work of Tesanovic, who employed 
a slave boson approach |36|. More recently, other studies of the effect of disorder in 
these model systems have appeared 

In order to make progress, we have applied the dynamical mean field theory of 
correlations and disorder to the Hamiltonian of Equation (3.1)[^ |o[ The 



derivation of the dynamical mean field theory is most easily accomplished on a 
Bethe lattice and we will focus on this class of models. In this case, the conduction 
electron density of states acquires a simple semicircular form. The solution of the full 
lattice problem then reduces to the solution of an ensemble of impurity problems 
supplemented by a self-consistency condition on the conduction electron Green's 
function ]4l|]. More precisely, the impurity problem action for a random site j is given, 
on the Matsubara frequency axis, by 

^imp ^rj^Y^ [fliiLo^) + Ef + V^2A(i..„)) /^,(*c.„)] , (3.2) 
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where the infinite-U constraint is implied and 

A(*^n) = -— ^ ■ (3.3) 

Here, t is the conduction electron hopping matrix element and Gc(ti^) is the disorder- 
averaged local conduction electron Green's function. Note that A(ia;„) is the 
hybridization function of the conduction electrons which is "seen" by each local f- 
site. Therefore, it corresponds to the local conduction electron Green's function with 
one f-site removed. This is the so-called "cavity" Green's functionpl|. 
The self-consistency condition determines G'c(w) through 



where 



GS^n) = {—- — -— — , (3.4) 



^A^^n)^- ^/^'^i^p.. (3.5) 



Here, (. . .)^^ denotes the average over disorder defined by the distribution functions 
Pi{Ef) and P2{V). S^P(itt'„) is the f-electron self-energy derived from the impurity 
model of ( [3.2[ ). This is explicitly defined by 



S}7(jc.„) = zLUr. - Ei- - Vf/^iiu^) - G}7(iu;„), (3.6) 
where the f-electron Green's function is given by 
G^7(r) =-(T/V(r)/t(0))7P 

G/7'(*^«)= /''rfrG;7(T)e^-"^ (3.7) 



Here, we have used (. . .) ''"^ to denote the quantum mechanical/thermal average under 
the action of Equation ^3.21 ) . 

Once the problem defined by (3.2)-(3.7) has been solved, the conduction electron 
self-energy Ec(*'^n) is obtained from 



2t 



GciiiOn)^ de . '''r V Po(e) = -v/l-(£W- (3.8) 

This self-energy is the important object for the calculation of the conductivity, which, 
in the infinite coordination limit, involves no vertex corrections p2| 

.^jpLr r /c+"),„(,)^(, , ^ (3.8, 

where /{v) is the Fermi function and we have introduced a lattice p aram eter a and 
the correct dimensional factors appropriate for three dimensions. In (3.9), A{e,u!) is 
the conduction electron one-particle spectral density 

A{e, u;) = ImG,(e, co) = Im (— ^ ^ . . (3.10) 
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The DC conductivity is then given by 



This expression can be further simphfied as shown in Appendix A . To make contact 
with the usual Drude formula for the DC conductivity it is useful to define the 
scattering time 

r ^ dv (-^\ depo{e)AHe, v), (3.12) 

which reduces to the usual case when t ^ 1/D 



In Appendix B, a formula which expresses the conduction electron self-energy in terms 
of the disorder-averaged impurity T-matrix is derived. This formula will be useful in 
Section ^, when we analyze the temperature dependence of the resistivity. 

Thermodynamic properties within this approach are not drastically different from 
the predictions of the simple disorder model of Ref. discussed in Section ||. 

Though the self-consistently determined conduction electron density of states "seen" 
by the f-sites (A(a;)) can be substantially modified as compared to its free value, these 
modifications only lead to a renormalized density of states, which is the quantity that 
ultimately enters into the Kondo temperature expression. Once the renormalized 
density of states is given, the prediction of the dynamical mean field theory is 
essentially the same as that of the simple disorder model of Ref. the difference 
being negligible. We discuss the thermodynamics properties in the dynamical mean 



field theory in Appendix C 



We would like to emphasize at this point what processes are included and what 
processes are left out of this approach. When the interaction U is turned off, the 
treatment of the disorder problem that is obtained is equivalent to the well-known 
coherent potential approximation (CPA)|^3|. This approximation is known to give 
reliable results as long as localization effects are negligible, since the conduction 
electron density of states fluctuations are treated on the average. This seems to be a 
safe approximation in the case of the alloys of interest, where estimates based on the 
zero temperature DC resistivity give kpl ~ 3 — 10. However, the interplay of disorder 
fluctuations and correlations on the f-sites is fully kept in our treatment and, indeed, 
is at heart of the physics that will emerge. This is evidenced by the fact that one 
needs to correctly solve an ensemble of interacting impurity problems in order to close 
the equations. Each member of the ensemble of impurity problems has an associated 
characteristic Kondo temperature Tk, thus generating a distribution of local Kondo 
scales H^. The fluctuations associated with this distribution of Kondo temperatures 
will be responsible for the anomalous low-temperature behavior in the non-Fermi 
liquid regime. Other processes which are left out of this essentially local approach 
are related to the RKKY interaction between f-sites mediated by the exchange of 
conduction electron spin fluctuations. We will comment on the possible limitations of 
this approximation in Section |^. 
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4. The zero temperature state 

The most difficult part of solving ([3.2|)-(3.7) is solving the impurity problem. Let us 
consider initially the clean case. Then, there is only one impurity model to be solved. 



Equations (3.4) and (3.8) are then trivially solved and yield 'Sc{uj) — ^{uj). At low 
temperatures, the impurity problem is governed by a Fermi liquid fixed point ^ and 
S'i"''(tj) can be parameterized at low temperatures and energies as [^^ 



T.fP{u;)^a + buj + ic{uj^ + n^T^) , (4.1) 
where, a, b and c are constants. From this, it follows that 

Sc(a;)« z 2^ (4.2) 

UJ ~ ef — izc (w^ + TT^i ^) 

where we have redefined a and b in terms of a renormalized f- level energy e/ and a 
wave function renormalization factor z. The first thing to notice in this expression is 
that at T = 0, Sc(0) is real, leading to a purely real conduction electron self-energy 
in the pure case. This is the hallmark of the low temperature coherent transport of 
the clean system, a consequence of translational invariance and a feature naturally 
incorporated in the dynamical mean field theory. Furthermore, the Fermi liquid form 
of the impurity problem self-energy ultimately leads to the Fermi liquid behavior of 
the conduction electron lattice self-energy. In particular, the resistivity derived from 



(4.2) will exhibit a characteristic law. 

Let us now move on to the disordered case at zero temperature. We have employed 
the slave boson mean field theory to solve the ensemble of impurity problems at 
T = 0^^. This theory is known to provide a good description of the infinite-U 
Anderson impurity model at temperatures T <ti Tx- It will also serve as a good 
starting point for the understanding of the complete solution of our dynamical mean 
field equations at an arbitrary temperature. Details of the slave boson treatment of 



the impurity problem are given in Appendix D. These solutions to the ensemble of 



impurity problems were then used in an iteration scheme to solve the full set (3.2)- 
(3.7). 

Figure Q shows the self-consistent solution to G'c(w) in the clean case and Figure Q 
shows the corresponding "cavity" Green's function A(cl'). The only modification 
introduced by the lattice of f-sites on the conduction electron density of states is 
the appearance of a gap centered around the renormalized f- level position ef. This is 
familiar from the slave boson large-N solution of the infinite-U Anderson lattice jssl]. 
The scale of the gap is given by w r^V^/D, which is the Kondo temperature scale. 
The important feature of the imaginary part of the "cavity" Green's function A(w) is 
the appearance of a delta function at the center of the gap e / , which corresponds to 
the removal of one f-site. 

The presence of disorder introduces important changes in these structures at the 
scales given by the range of Kondo temperatures produced as a consequence of the 
distribution of local f-parameters. These are defined here as 



TK,^y'(eJ-A;.(o))2-i-(A;'(o))2, (4.3) 

where Aj(a)) = r'jV^ A{uj) and we are denoting real and imaginary parts by single and 
double primes, respectively. In the Kondo limit, they acquire the familiar form 



T, 



De^^(Ef/{2poVf)), (4.4) 
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Figure 4. (a) Imagi nar y and (b) real part of the conduction electron "cavity" Green's 
function A{u>) (see ( |3.3| )) as a function of frequency at T = 0. The parameters used 
were: D = 1, V = 0.1, fj. = 0, = —0.05 and there is no disorder. 
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where po = A"(0)/7r. Figure ^ shows a typical distribution of Kondo temperatures 
from a fully self-consistent solution of the zero temperature problem. For comparison, 
we have also plotted the P{Tk) used in the fit to the susceptibility of UCu4Pd 
in Ref. pO[ |. The self-consistent distribution has a structure very similar to the 
experimentally determined one. 

60 I ' 1 ' 1 




' ' ' ' ' 

0.00 0.02 0.04 

Tk/d 

Figure 5. Typical distribution of Kondo temperatures obtained in the fully self- 
consistent solution of the dynamical mean field theory ( ) and the distribution 

appropriate for UCu4Pd from Ref. |3o[ ( ). The distribution was a Gaussian and 

the parameters used were: D = I, n = = -I, {V'^) = 0.17, Wy2 = 0.032. The 

upturn at very low Tk's is not shown. 

These modifications are illustrated in figures ^, ^ and which show the influence 
of various amounts of disorder in the gap structure of Gc(w) and A(cl;). Only the low 
energy region is shown since that is the only region that is strongly influenced by the 
introduction of disorder. It is clear from Figure ^ that sufficient disorder can close the 
gap in the conduction electron density of states. Qualitatively, a distribution of Kondo 
energy scales for the ensemble of impurity problems will be reflected in a distribution 
of energy gaps. Their sizes and positions are essentially determined by the Txj^s and 
e^'s, respectively, and sufficient disorder in their distributions will lead to an ultimate 
closing of the gap. In the metallic regime being considered here, the disappearance 
of the gap does not lead to dramatic effects, since the chemical potential lies away 
from the gap region. However, there might be important consequences in the Kondo 
insulator or Anderson insulator regimes. Corresponding modifications are also present 
in the real parts of Gc{uj) and A(w), as seen in figures ^ and ||. The result is similar 
in the case where is uniform and Vj is distributed, since the relevant scales are the 



15 



iKj on E'^ 



f 



Kondo temperatures (4.4). However, due to the different dependences of Tj 
and Vj, their distributions will be quantitatively different and will reflect differently 
on the Tr-j distribution and consequently on the transport properties. 
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Figure 6. Imaginary part of the conduction electron Green's function as a function 
of frequency at T = for different amounts of disorder in the E-^ parameter. A 

uniform distribution was used with widths W = 0.001 ( ), W = 0.05 ( ) 

and W = 0.15 ( ). The parameters used were: D = I, V = 0.1, /i = 0, 

(Ef) = -0.05. 



The central aspect of the interplay between disorder and correlations is an 
enhancement of the bare disorder by the local Kondo physics of the impurity ensemble. 
This can be quite easily understood by an examination of the qualitative aspects of 
the solution at T = 0. At zero temperature, only $j(0) enters into the calculation 
of Sc(0) or Gc(0), either of wh ich, in turn, is enough for the calculation of the 
DC conductivity (see (|]T|) or (^)). As was emphasized before, given a certain 
distribution of the quantities $j. Equation (3.4) corresponds to a CPA treatment of 



disorder in the conduction band, where the <i>j's play the role of scattering potential 
strengths. Interactions are important in determining the value these latter quantities 
acquire, given a self-averaged conduction electron bath through Aiuj). 

From the Fermi liquid analysis, the zero temperature, zero energy form of $j(0) 
is real and given by 

$(0) = , . (4.5) 




Figure 8. Real part of the conduction electron "cavity" Green's function A(uj) as a 
function of frequency at T = for different amounts of disorder in tlie parameter. 

A uniform distribution was used with widths W = 0.001 ( ), W = 0.05 ( ) 

and W = 0.15 ( ). The same parameters were used as in Figure H. 
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In the slave boson mean field theory, it is given by (see (D.ll)) 

„2t/2 



J ^3 



(4.6) 



Now, in the Kondo limit >> PoV^ {Ej < 0), it is easy to show from the mean 
field equations (D.7) and (D.8) that 



A'(0) 



(4.7) 




0.05 



0.15 



0.25 



-0.15 -0.10 -0.05 0.00 



Figure 9. Variation of the effective local scattering potential strength 3>(0) as a 
function of and . Notice how <I'(0) varies on the scale of conduction electron 
parameters (D) not on the scale of local f-parameters. The distributions used were 
Gaussian and the parameters were D = 1, = 0, (a) E-f = —1, (V^) = 0.17, 
Wy2 = 0.032, (b) V = 0.1, (Ef) = -0.055 and Wj^f = 0.027. 

The important point here is the fact that the impurity problem parameters have 
disappeared from the local effective scattering potential strength and its scale is now 
given by the conduction electron band scale. As can be seen from Figure ||, the 
real part of the "cavity" Green's function that enters (4.7) spans a wide range of 
values in the region close to the chemical potential (here set to zero). Because of 
the distribution of local parameters and rj, this whole region close to /x will be 
probed by the diflFerent f-sites and <&j(0) will also vary. However, its variation will be 
on the conduction electron scale. Indeed, we have plotted in Figure the variation 
of ^j{0) for a given strength of disorder, in both cases of random Eps and random 
Vj 's. It is clear that its distribution range is not on the same scale as the variation of 
local f-parameters, which is very narrow, but rather, on the scale of D. Therefore, the 
effective disorder seen by the conduction electrons is considerably enhanced due to the 
local f-shell correlation effects. 

As a result of this enhancement effect, large resistivities can be generated even 
though the range of variation of the f-shell parameters is very narrow. This is 
illustrated in Figure |l^, where the scattering rate l/r, as given from ( 3.12| ), is 
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Figure 10. Scattering rate as a function of the width of the uniform E-^ distribution. 
The strong correlations in the f-shell produce an enhanced effective disorder. The 
parameters used were: D = 1, fj, = 0, V = 0.1, (E^) = —0.05. 



plotted as a function of the width of the uniform distribution. It is clear that 
modest amounts of disorder in i?-^ produce rather substantial scattering rates. Thus, 
although the clean system exhibits coherent transport at low T, the introduction of 
a small amount of f-site disorder leads to the destruction of these coherence effects. 
Qualitatively, once the lattice effects of coherence are destroyed by sufficient disorder, 
the resistivity as a function of temperature will then resemble the independent 
Kondo impurity results, with its characteristic decreasing resistivity with increasing 
temperature. This point was nicely illustrated in a recent study ||3^ of binary alloy 
disorder in the Anderson lattice, which relies on a similar treatment of correlations 
and disorder. Furthermore, several doping studies of heavy fermion systems seem to 
bear out the above picture [p^. 

We summarize now the two major results of the study of disorder discussed in this 
section: 

• Due to the local Kondo physics at each f-site, the effective disorder generated 
from a bare distribution of local f-shell parameters is strongly renormalized up to 
scales of the order of the conduction electron bandwidth. 

• Although the clean system has low resistivities due to the onset of coherence at 
low T, moderate amounts of f-shell disorder are capable of destroying this low-T 
coherence, leading to characteristic incoherent Kondo scattering behavior. 
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5. The linear temperature dependence of the resistivity 

Having established the incoherent nature of the transport with sufficient disorder 
strength, we now focus on the temperature dependence of the resistivity in this 



strongly correlated disordered state. We will rely on Appendix B, which relates 



the conduction electron self-energy T.c(lu) to the averaged impurity model T- matrix 
T"^^{u!) (Equation (B.4)). At sufficiently high temperatures, compared to the highest 
Tk in the distribution, the averaged T-matrix at low energies becomes very small, 
reflecting the weak coupling nature of the impurity response at high temperatures. 
Therefore, there will be only a small contribution to the conduction electron self-energy 
and consequently to the resistivity. At zero temperature, as discussed in Section ^, 
sufficient disorder will kill all coherence in the transport properties and the system 
will exhibit a rather large resistivity. What happens at low temperatures? 

For that, it is convenient to analyze the average impurity T-matrix. We have 
therefore plotted in Figure the imaginary part of the impurity T-matrix, using 
the distribution of Kondo temperatures derived from the fit to the thermodynamic 
properties of UCu5_2;Pda;[0- It is clear that the leading temperature dependence is 
linear. This leads to a linear conduction electron self-energy and to a linear resistivity 



through Equation (3.11), consistent with the behavior observed in several of the alloys 
in Table 0. At higher temperatures, there are clear deviations from the leading linear 
behavior as can be seen in the inset. Can we understand what conditions are necessary 
for this anomalous non- Fermi liquid behavior? 

The essential condition for the linear temperature dependence of the resistivity, 
like the thermodynamic response discussed in Section |2[ is that the distribution of 
Kondo temperatures be such that P{Tk = 0) ^ 0. In other words, there must be 
a finite fraction of f-sites with arbitrarily small Kondo temperatures. This can be 
made more apparent by an analysis similar to the one in Section |^. However, one has 
to be careful of how one understands the averaging process. This is because, unlike 
the thermodynamic responses, transport properties cannot he averaged over single- 
impurity results. This is made evident in the analysis of the zero temperature state of 
Section ^. The single impurity result amounts to a strong, almost unitary scattering 
center at T = 0, with a correspondingly large resistivity. Any attempt to simply 
average this result would lead to a large resistivity. However, due to translational 
invariance, it is clear that the clean sample can only have a vanishing resistivity, a 
result which encounters a natural description in the dynamical mean field theory, 
as was explained. Even in the disordered case, only the CPA-like averaging process 



embodied in Equation (3.4) has any sense. However, given a disordered ground state 
deviations from zero temperature can still be analyzed in a fashion that resembles the 
averaging of thermodynamic quantities. This is because, at low temperatures, a few 
low-T/f spins are unquenched and cease to contribute significantly to the scattering. 
Since they form a dilute system of subtracted scattering centers, their contribution is 
additive and can therefore be averaged. 

With these caveats in mind, we can proceed to analyze the conduction electron 
self-energy. If the temperature is raised from to T, there will be corresponding 
variations in all quantities. From (B.3) 



l-t^Gliu:) 

d2.c(t^) = 



5Gc{io) (5.1) 

=0 
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Figure 11. Temperature dependence of the imaginary part of tlie single impurity 
T-matrix averaged over the disorder distribution appropriate for UCU3 sPdi.s, as 
determined experimentally in Ref. | po| . The inset shows the same quantity over a 
wider temperature range. 



and SGc{uj) can be obtained from (B.2), 

{^2 _^ _ t^Gciu;)] [iJ - 3t'Gc(a;)] SGdu;) = 5(r;'"P(a;))^\ (5.2) 



At t his point it is useful to remember the definition of the T-matrix (see (BJ^) and 

(H)) 



(5.3) 



In (|5.3| ) we have highhghted the fact that is a functional of Gc(ti-') as well as 

a function of w and T. Therefore, it is clear that the T-matrix depends on the 
temperature through Sy™^. However, since Gc(ti-') changes with temperature, we also 
need to take into account this additional contribution, both explicitly and implicitly 
through '^Yj' ■ therefore, separate three different contributions 



5r;-p(u;) = &{rf^{uj) + JsTf p(w) -f hTf'ioj), 



(TfP(c.))" 



<^irr'(^) 



<5Gc(c^), 



T=0 



(5.4) 
(5.5) 

(5.6) 



21 



y2 

3 



dto' 



(5.7) 



T=0 



In (5.6), we have used St to denote the variation with T which is not imphcit through 
Gc{co) and G° is used as a reminder that the variation is calculated with a fixed 
zero-temperature bath. Equation (5.7) involves the functional derivative of the f-self- 
energy with resp ect to the bath Green's function. Equations (5.4)^(5.7) can now be 
substituted into (5^) to yield 



t'^+[uj~ t'^Gciuj)] [to - St^Gcito)] 



i2( r-P(a;) ) 



SG.iLo) 



T=0 



SGJuj 



'lo - t^Gc{Lu)Y 

SG,W) = (<5tT;"P(^))^ 



GO 



(5.8) 



T=0 



This is an integral equation for 5Gc{lo) whose source term is given by the averaged 
variation of the T-matrix with the temperature. Without the last term on its left-hand 
side, Equation (5.8) is actually a simple algebraic equation. The last term describes 
the feedback effect that the change in the self-consistent conduction electron bath 
generates on the ensemble of local impurity actions. Since raising the temperature by 
a small amount leads to the unquenching of a few dilute spins, we do not expect this 
feedback effect to be large and will thus neglect the last term on the left-hand side of 
(5.8). 

We now analyze the source term on the right-hand side of (5.8). At zero frequency, 
the imaginary part of the T-matrix can be written as 



Im7;™P(T) 



sin^ i5i 
T^Po 



(5.9) 



where 6oj is the phase shift at T = 0. We can write down the asymptotic limits of the 
scaling function t{x) 



t{x) 
t{x) 



1 



- ax 



(a;< 1); 
(a;> 1), 



(ln(2;))2 

where a and /? are universal numbers. Therefore, 



5Timr;'°p(T) 



sin^ Soj 



1 



T 



sin^ (5( 



(5.10) 



F{T/Tk,1 (5.11) 



which defines the function F(x). This function will be averaged over with the 
distribution of Kondo temperatures. Therefore, keeping T fixed and as as function 
ofTK 



F{T/Tk,) 



F{T/Tk,) « 1 



Mt/Tk,W 



{Tkj » T); 



(5.12) 
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It can be seen that F(T/Txj) has a peak at Tfcj — with width T, decaying rapidly 
to zero as 1 /T^^ for large Txj (Figure |l^) . For low temperatures compared with the 
typical scale of the distribution function P{Tk) one can write 



I sin 6( 



'Oj 



imp 



(5.13) 



where a — J dxF{l/x). A similar analysis can be carried out for the real part of the 
T-matrix. Therefore, after averaging over Tx one gets 



(<5TlmT;'"P(^)) 



aPn sin Sf) 



-T, 



(5.14) 



consistent with Figure |T^. It is clear that, as long as the distribution of Kondo 
temperatures has finite weight at Tr- = 0, the average T-matrix will show a linear 
temperature dependence. If P{0) = or negligible, then Fermi liquid behavior is 
recovered, with the characteristic law. The result of (5.14) should be plugged into 
(5.8) and then into (pM for the final expression of the conduction electron self-energy. 
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Figure 12. Experimentally determined distribution of Kondo temperatures of the 

alloys UCu5_^Pdi, with x = 1( ) and 1.5( ) (from Ref. and the 

function F{Tk,T) ( ) defined in the text. The upturn at very low 'i'j^'s is not 

shown. The function F{Tk, T) only probes the Tk = value of the distributions at 
low T. 



As our analysis above has shown, the same physical mechanism underlies 
the anomalous non-Fermi liquid behavior of the resistivity as in the case of the 
thermodynamic properties. It is the presence of 1ow-Tr- spins, unquenched even 
at low temperatures, which gives rise to anomalous scattering. Though the zero 
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temperature transport is a reflection of the full structure of the distribution function, 
the leading low temperature behavior is much simpler, corresponding to the gradual 
decrease of the number of effective scattering centers. Since the number of released 
spins is small at low temperatures, their effect is additive and an average over their 
subtracted T-matrices is well justified. Finally, we add that an immediate consequence 
of the physical origin of the anomalous scattering in this disorder model is a negative 
magnetoresistance at low temperatures. Much like the temperature, a magnetic field 
acts to destroy the low-Tfj- Kondo singlets and thus to suppress their effectiveness as 
sources of disorder. 



6. The dynamic susceptibility 

The dynamic susceptibility of VCu^-xPdx (x = 1 and x — 1.5) has been measured 
with inelastic neutron scattering and reported in Ref. [ p3| . The first important 
result was that the (f-dependence of the magnetic response could be completely 
accounted for by the g-dependence of the Uranium ion form factor, suggesting that 
the spin dynamics is completely local. Furthermore, as expected, the frequency and 
temperature dependences of the imaginary part of the dynamic susceptibility are 
anomalous. There was no significant difference between the magnetic behaviors of 
the X = I and the x — 1.5 alloys. These results provide a useful testing ground 
for the disorder model. Indeed, one can use the distribution of Kondo temperatures 
deduced from the fits to the thermodynamic data to determine the dynamic response 
and compare with the experimental data. Like the static magnetic susceptibility, the 
dynamic response will be dominated by a few unquenched spins. Therefore, it is a 
reasonable assumption to take the overall dynamic lattice response to be essentially 
given by an average over the single-impurity results. Moreover, the local nature of the 
measured dynamic susceptibility is consistent with this assumption. 

There is currently no complete description of the dynamic susceptibility of a single 
Kondo impurity for the full range of temperatures and frequencies, though the methods 
to carry out this task certainly are available. Among the existing results, we cite the 
unpublished work of Costi and Hewson for the dynamic susceptibility of the Anderson 
model, quoted in Ref. |Q. Besides, the non-crossing approximation (NCA) and an 
extension of it have also been used to determine this response p^. Quite often, the 
dynamic susceptibility of a Kondo impurity is fitted to a relaxational formp7[| 

f,, xiT)T{T)uj 

^ ("'^^^ c.^ + r^(r) ' (^-'^ 

where x(^) is the impurity spin susceptibility and the linewidth r(T) is a function 
of temperature only. Note that this form automatically satisfies the Kramers-Kronig 
relation 

TT J-oo ^ 

The high temperature behavior of r(r) is given the Korringa law with logarithmic 
corrections 

r(T) » 47r (poJf T [1 - 4 (poJ) InT] . (6.3) 
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At zero temperature, one can use the so-called Shiba relation to determine r(0)p^ 

lim^>^^^=.r(0)^^ (6.4) 

where w « 0.4107 is the so-called Wilson number. In the absence of a better 
description, we have employed a crude approximation that interpolates between these 
two limits 

r(T)=<; ^ ^ for T<Tk ^ 

47r(poJ) T + a for T > Tk 

where a is such that r(r) is continuous. We stress the crudeness of the approximation 
and regard it as a rough description of the actual behavior. 
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Figure 13j_^ Comparison between the measured dynamic magnetic susceptibility 
from Ref. p3 and the prediction of the disorder model. Here, we have used the 
distribution function appropriate for a; = 3.5 according to Ref. pcl]. 

Once the behavior of the dynamic susceptibility for a Kondo impurity is assumed to 



be the one given by Equations (|6.l|) and (6.5), one can then easily perform the average 



with the distribution of coupling constants determined experimentally in Ref. | |30| . We 
have done so and the results are shown in Figure |l| (a; = 3.5) and Figure |l| (a; = 4). 
We stress that, once the distribution off-shell parameters is determined from the fits to 
the static magnetic susceptibility, no additional fitting is performed. As can be seen 
from the figures, the agreement between the experiment and the predictions of the 
disorder model is rather good, considering the range of frequencies and temperatures 
and the crudeness of our assumptions. The agreement is slightly better when one 
uses the distribution of the x = 3.5 alloy. Though a more accurate description should 



25 



be endeavored, we believe the disorder model cannot be ruled out by the neutron 
scattering data. 
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Figure 14. Comparison between tlie measured dynamic magnetic susceptibility 
from Ref. and tlie prediction of the disorder model. Hem, we have used the 
distribution function appropriate for x = 4 according to Ref. |30| . 



7. Discussion and conclusions 

We would like to now pause and consider the overall picture that emerges from the 
disorder model as well as the drawbacks of our current treatment of the problem. 
We have emphasized throughout that the important ingredient at the base of all the 
non-Fermi liquid features in this approach is the presence of f-sites with arbitrarily 
low local Kondo temperatures. More precisely, P{Tk = 0) 7^ 0. As the temperature is 
raised, these low-Tjf spins are gradually unquenched and it is this very release process 
that gives rise to the anomalous thermodynamic as well as transport properties. 

One of the practical difficulties one faces in applying the dynamical field theory, 
is the fact that, due to the very nature of the physics involved, one needs to be able 
to solve the ensemble of impurity models in the whole range of temperatures from 
T ^ Tkj to T ^ Tkj ■ This is a notoriously difficult task as is evidenced by the long 
history that led to the final solution of the Kondo problem. The conventional approach 
of Wilson's numerical renormalization group, as is currently formulated, relies on a 
special "energy-shell" decimation procedure, which is not obviously adapted to the 
present case of a general conduction electron density of states. Therefore, in the non- 
trivial case of the resistivity, we are not yet able to reliably predict the value of the 
coefficient of the linear term, which might be checked against experiments. 
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The dynamical mean field theory does not take into account the RKKY interaction 
between f-sites. These have been studied most extensively in the context of the two- 
impurity Kondo problem^^. From these studies it is known that, if the RKKY scale 
is large enough compared to Tk, the two impurity spins develop strong correlations 
which tend to lock them into a singlet state, whereby the Kondo effect is killed. Since 
we rely on the presence of certain sites with very low Kondo temperatures, one might 
wonder whether the inclusion of RKKY interactions would not effectively provide a 
low energy cutoff below which P{Tk) would be essentially negligible. Although we 
cannot give a rigorous answer to this question, some very general arguments indicate 
that the possibility of having 1ow-Tr- spins is more robust than one might expect. 

We have already emphasized that, if the distribution function P{Tx) is broad 
enough, at low temperatures compared to its overall width, the fraction of spins 
which remain unquenched is rather small. Therefore, in general, they form a dilute 
system of spins of density niow , the average distance between them being proportional 
to {niow)~^^^ ■ Since the average distance will be large, this will render the RKKY 
interaction less effective, since its strength decays as R~^. Furthermore, we will argue 
that the effectiveness of the RKKY interaction in suppressing the Kondo effect is 
actually a higher order effect in the dilution niow 

Indeed, consider two Iow-Tr- spins chosen at random in the sample. In general, as 
argued above, they will be far apart. Now, unlike the usual two-impurity Kondo 
problem, these two spins will have different coupling constants Ji and J2 and, 
therefore, two different Kondo scales Tki and Tk2- The RKKY scale Trkky 
J1J2/D, apart from the dependence on the distance between the spins. Now, in the 
conventional two- impurity Kondo problem, where Tki = Tk2i one has, qualitatively, 
two independent Kondo screening processes, one at each spin, when Tki — Tk2 ^ 
Trkky, whereas the effect will be suppressed when Tki = Tk2 Trkky- In 
our disordered case, RKKY dominates if, say, Tki < Tk2 ^ Trkky- However, in 
the intermediate case when Tki ^ Trkky ^ Tk2, two independent Kondo effects 
will survive, even though Trkky is still larger than one of the Kondo temperatures. 
Indeed, as one comes down in energy scale, by the time one hits Trkky i spin 2 has 
already undergone quenching and is no longer available to correlate into a singlet-like 
composite with spin 1. Therefore, the RKKY effectiveness depends on the random 
selection of two spins, both of which must have low enough Kondo temperatures. We 
would thus expect it to be of order rif^^- Naturally, the transitions between different 
regimes are all crossovers (in the absence of special unphysical symmetries) and there 
will be no sharp distinctions between the different cases and hence no sharp cutoff to 
the distribution function P{Tk)- 

We have emphasized throughout our analysis the importance of the structure 
of the distribution of Kondo temperatures, in particular, whether it intercepts the 
vertical axis at Tk = or not. Due to the exponential dependence of Tk on the 
f-parameters, very small modifications in the width of the bare distribution lead to 
rather large changes in the distribution of Kondo temperatures. If the bare width is 
too small, P{Tk) will be negligible at low Tk, whereas, if it is too large, P{Tk) will 
be divergent. Therefore, the situation where P(Tk) can be taken as approximately 
constant when Tk — > only holds in a narrow window of bare disorder widths, in our 
current approach. In a more speculative vein, one is tempted to think that the RKKY 
interaction, which is left out of the current approach, might intervene to impose an 
upper bound on the density of Iow-Tk spins which are allowed to exist in such a 
state. It would do so by effectively eliminating them through the formation of random 
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singlets 13 . 

Finally, by treating the disorder seen by the conduction electrons in a mean- 
field CPA-like fashion, we have neglected the effects of fluctuations in the conduction 
electron density of states. It is natural to expect that these eff'ects might act to further 
renormalize the distribution of Kondo temperatures. It would be interesting to study 
whether, like in many other treatments of disordered systems [p7[, these effects lead to 
a flow towards a universal form of the distribution function. 

In conclusion, we have presented a complete picture of the possible origin of the 
non-Fermi liquid behavior in Kondo alloys. Our study has shown the important effect 
correlations have on the role of disorder in f-electron systems, considerably enhancing 
the bare f-shell disorder strength. The interplay between disorder and correlations 
leads to the idea of a distribution of Kondo temperatures, whose structure determines 
whether a Fermi liquid description is possible. In the particular case when P(0) ^ 0, 
both thermodynamic and transport properties are anomalous and incompatible with 
a Fermi liquid picture. 
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Appendix A. Simplified expression for the conductivity 



In the case of a semicircular density of states, a simplification of the expression (3.11) 
can be achieved. Indeed [08|, 



4^Re 



dGi{z) 



dz 



ImG;(z) 
Im^; 



where z = w + /i — Sc('j-') and 

.Po(e) 



Gi[z) = J de^ 



Using the identity z = (z) + 1 Gi{z), which follows from (|3.8|), we get 



2e^ 



Re 



1 



{t^Gi{z)r-l 



1 



I ~ \f'Gi{z)\^ 



(A.l) 



(A.2) 



(A.3) 



Appendix B. Conduction electron self-energy and the impurity model 
T- matrix 

The expression for the conduction electron self-energy T,c{(jj) can be recast in a more 
illuminating form by employing the T-matrix associated with the impurity problem, 
which is defined by 

7j'"P(a;) = 1^'G™P(^). (B.l) 
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Inserting this definition into (3.6), then into (3.5) and ( p.4| ), we can write after some 
manipulations 



1 



(T™P(u;))-- 



Now, (|3.8| ) can be explicitly solved for a semicircular density of states giving 

Y.c{uj) =uj + ^ - — \ ^ - t^Gciio). 
Gc[oj) 

Combining (B.2) and (B.3) we finally get 



Sc(w) = = 



(T™P(c^))-- 



G,{uJ){L0 + ^i-t^G,{L0)) 



(B.2) 
(B.3) 

(B.4) 



Appendix C. Thermodynamic properties 

In this appendix, we will investigate whether the assumptions of the simple disorder 
model of Ref. |3^ (hereafter, SDM) are justified within the framework of the 
dynamical mean field theory. Let us recapitulate what the assumptions of the SDM 
are. Essentially, the model assumes a collection of independent Anderson or Kondo 
impurities, each of which hybridizes with the same, featureless conduction bath. 
The parameters of this collection of impurities are distributed according to some 
assumed distribution function(s). In the particular case of Ref. which assumed 
Kondo spins, the dimensionless Kondo coupling constant A = pqJ was assumed to be 
distributed according to a Gaussian. Thermodynamic properties are then calculated 
by taking the average over the single-impurity results with the appropriate distribution 
functions, as if the impurity responses were completely uncorrelated. 

It is clear that the dynamical mean field theory retains much of the flavor of the 
SDM, with the ensemble of impurity problems playing the role of the uncorrelated 
collection of spins. The question we would like to answer is whether the equivalence 
can be shown within the framework of the dynamical mean field theory and, more 
importantly, what are the conditions for the validity of the equivalence. We will 
confine our analysis to the total energy, which is enough for the calculation of specific 
heat. The same conclusions apply to the susceptibility, which, however, requires 
consideration of the free energy. 

Let us first write down the total energy in the framework of the dynamical mean 
field theory. We have 

ijcr ja 
JO- j 

We can write 

(4so> = 7^EGci(^^«)- (C.2) 

Now, when the coordination number z goes to infinity and for the case of a Bethe 
lattice (in the absence of symmetry breaking) one can prove that [p| 

G^i(zc.„) = -U,Gli^u;„), (C.3) 
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from which it foUows 



(C.4) 

where A/" is the number of lattice sites and, consistent with the correct rescahng of the 
hopping in the infinite coordination hmit, = ztfj. Plugging (C.4) into (C.l), one 
then has the total energy completely expressed in terms of local quantities. This is 
the quantity we want to compare with the prediction of the SDM. 

In order to do that, one needs to define what precisely one means by the SDM. In 
the dynamical mean field theory, we have the ensemble of impurity models defined by 
Equation (3.2). We, therefore, define the total energy in the SDM by the sum of the 
energies of the various impurities 



(C.5) 



To define the energy of one impurity, we first write the Hamiltonian corresponding to 
Equation (O) 



j 



ki7 



k(7 



where we introduced fictitious fermionic operators a^^ to mimic the self-consistent 
conduction electron bath through the "cavity" Green's function of Equation ( |3.3| ) 



A(c.) ^ E — 



We then define the appropriate single-impurity energy by 



(C.7) 



(C.8) 



kij 



where we subtracted the total energy of the fictitious conduction electron bath 
evaluated at y = 0. Now, the fictitious conduction electron Green's function is 



Using 



Ga{k,k',u:)^5j:-^,Gl{k,uj) + Gl{k,uj)V^Gf{iu)Gl{k\uj). (C.9) 
= TY,Ga{k,k,iujn), (C.IO) 
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we have 
Thus 
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= Tl^2^G^(io.„) \Gl{k,iuj^) 
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In the last equality we used Equation (p.7| ). 

We can now use Equations (C.l), ( p.SD , ( p.6D and (C.8) to write the difference in 
energy between the SDM and the dynamical mean field theory 



AE = E 



DMF 



-E. 



SDM — 



n 



k(T 



Et 



^ ka ka' 



— [a- a- 

y— Q ka ka 



To simplify things further, we can use Equations (p.4D and (C.12) to get 



AE = 2Tj2{^t^Gl{i 



A(za;„) + iio, 



dA{iujn) 



d{iuj„) 



We then use Equation (B.2) and the definition ( |3.3| ) to arrive at our final expression 
for the difference in energy between the dynamical mean field theory and the SDM 



(C.14) 



AE 



Gc{iuJn) 



A{iu}n) 



- 1 



1 



d\nA{iojn) 



d\n{iujn) 



(C.15) 



We are now in a position to determine whether the SDM is accurate enough to give 
the thermodynamic properties of the system within the dynamical mean field theory 
framework. For this, it is enough to consider Equation (C.13). All the quantities in 
this equation are related to conduction electron kinetic properties. The corresponding 
densities of states are given by the imaginary parts of either Gc(w) or the "cavity" 
Green's function A(w). As can be clearly seen from Figures ||and^, these quantities 
are weakly renormalized. Their contributions to, say, the specific heat coefficient 
7 = Cv{T)/T are of the order of 1/Z? and are completely negligible when compared 
to the contribution from the impurity part, which is of order 1/Tk 3> 1/ D. Therefore, 
as far as thermodynamic properties such as 7 and x are concerned, the approximation 
of averaging over single-impurity results is perfectly consistent with the solution of 
the full dynamical mean field theory. This is true of both the clean and the dirty 
systems, since nothing in this argument relied on the presence of disorder. However, 
we stress the fact that the ensemble of impurity problems must be solved in the fully 
self- consistent conduction electron bath. Though the renormalizations of this bath are 
small, they affect the impurity properties through the conduction electron density 
of states po- The latter quantity appears in the argument of the exponential in the 
expression for the Kondo temperature and can, therefore, lead to rather substantial 
renormalizations of the total energy. To fully describe these changes, one needs to 
solve the full dynamical mean field theory equations. 



Appendix D. The slave boson mean field theory of the impurity models 

The slave boson description of the infinite-U Anderson Model starts with the 
replacement of the non-holonomic constraint Uf < 1 imposed by the [/ — > 00 condition 
by a new set of slave boson operators hj together with a holonomic constraint [35| 



fja 

n-fja 



n-fja 



b]b, = l. 



(D.l) 
(D.2) 
(D.3) 
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The constraint is then imposed by introducing a Lagrange multipHer term in the 
Lagrangian 



Co + i drXj (t) 



(D.4) 



where Aj(r) is an additional bosonic field variable. The transformed action 
corresponding to the impurity action (3.2) is then 



rtimp rrimp , oimp 

47= ["dr f''dr'Y,fl{T)b^iT)V^^A{T-r%^{r')b]{r% 
Jo Jo 



(D.5) 



where A(r) is the Matsubara-Fourier transform of ( p.3[ ). Note that the second part of 
the impurity action, which corresponds to hybridization processes between the f-site 
and the effective conduction electron bath, has been modified by the introduction of 
the slave boson field 6j(t), which makes the bookkeeping of the occupation of the 
f-site. 

At the mean field level, the bosonic fields 6j(t) and Aj(r) acquire a time- 
independent expectation value and behave as c-numbersjs^. The resulting effective 
action is quadratic in the pseudo f-electrons and can be exactly solved. Following the 



convention of writing i{Xj) — ^ E!j and (bj) = Vj, it reads 

Sf^Tj2 (-*^n + + r]V^^A{zu;^)) /^,(zc.„)] +(e^^ - i?/) (r^ - l) .(D.6) 

The mean field parameters rj and e - are determined by a saddle point extremization 



of the free energy corresponding to ( |D.6D . The mean field equations then read 

/(^)A;-(c.) 



du) 

DO 

duj 



J 



^-1 



0; (D.7) 



f{L0)A''{uj){L0 



{ef-Ef) =0, (D., 



where /(w) is the Fermi function, Aj{uj) — r'jV^A{uj) and single and double primes 
denote real and imaginary parts, respectively. We note that, given the hybridization 
function A{iujn), each value of the bare parameters Ej and Vj will define a different 
impurity problem with its own corresponding values of rj and Cj . Thus, a distribution 
of mean field parameters is also generated. 

The mean field treatment we have described gives the following expression for the 
f-electron Green's function 



Gfj{iuJn) 



i^n - ej - r]VfA{iuj^) 



(D.9) 



32 



Note that the numerator in (D.9) is crucial. It is a consequence of the slave boson 
prescription (D.l) and distinguishes the pseudo f-electron Green's function, from which 
it is absent, from the real f-electron Green's function. Finally, using the definition 
(3.6), one can write the mean field expression for the f-electron self-energy 



El 



(D.IO) 



which yields 



3 ^ ] 



(D.ll) 



This can then be inserted into (3.4), thus closing the self-consistency loop. Note that, 
in the pure case. 



I],(w) = $(c^) 



(D.12) 



which is consistent with the low-energy Fermi liquid parameterization of (4.2). 
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